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ABSTRACT 

The  oscillatory  forces  generated  by  the  interaction  of 
a  train  of  gravity  waves  with  a  rigid  horizontal  circular 
cylinder  located  near  a  simulated  ocean  floor  were  deter- 
mined experimentally.   A  circular  cylinder  of  six  inch  dia- 
meter was  subjected  to  gravity  waves  of  wave  lengths  varying 
from  2  to  18  feet  in  a  two-dimensional  wave  channel.   Both 
the  horizontal  and  vertical  components  of  the  hydrodynamic 
force  were  measured. 

Experimental  results  are  plotted  in  dimensionless 
form.   The  horizontal  and  vertical  dimensionless  force  co- 
.  efficients  are  plotted  as  a  function  of  the  relative  water 
u.eptn,  trie  relative  wave  length,  anu  the  relative  wave 
height.   It  was  found  that  the  horizontal  force  coefficient 
varied  quite  linearly  with  the  wave  height,  and  the  maximum 
value  was  approximately  equal  in  either  direction.   The 
vertical  force  component,  however,  showed  a  much  more  com- 
plex variation.   The  upward  and  downward  magnitudes  were 
approximately  equal  for  small  amplitude  waves,  but  as  the 
wave  height  increased,  the  variations  became  quite  non- 
linear and  the  magnitudes  in  the  two  directions  differed 
considerably . 
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NOMENCLATURE 

Symbol  Description 

a  radius  of  circular  cylinder 

C  added-mass  coefficient,  dimensionless 

d  depth  of  submergence  of  cylinder 

F  maximum  horizontal  force  per  unit  length 

f  F  /  [pga2 (H/2a) ] ,  horizontal  force  coefficient, 


x 


dimensionless 


F  maximum  positive  vertical  force  per  unit  length, 

y  (upward) 

f  F  / [pga2 (H/2a) ] ,  positive  vertical  force,  coeffi - 
y      y  cient,  dimensionless 

F  maximum  negative  vertical  force  per  unit  length, 

^  (downward) 

f  F  / [pga2 (H/2a) ] ,  negative  vertical  force  coeffi - 
^      ^  cient,  dimensionless 

g  acceleration  due  to  gravity 

H  height  of  the  incident  wave 

h  water  depth 

k  wave  number=  2tt/L 

L  wave  length  of  the  incident  wave 

1  length  of  the  cylinder 

T  period  of  the  incident  wave 

t  time 

u  component  of  fluid  velocity  in  the  x-direction 

u  component  of  fluid  acceleration  in  the  x-direction 

v  component  of  fluid  velocity  in  the  y-direction 

v  component  of  fluid  acceleration  in  the  y-direction 


Symbol  Description 

x  space  coordinate 

y  space  coordinate 

$  velocity  potential 

p  density  of  the  fluid 

a  frequency  of  the  incident  wave=  2tt/T 

y  dynamic  fluid  viscosity 


DIMENSIONLESS  NOTATION 


d/a  relative  depth  of  submergence 

F  / [pga2 (H/2a) ]  horizontal  wave  force  coefficient 

F  /  [pga2 (H/2a) ]  vertical  wave  force  coefficient 

H/2a  relative  wave  height 

h/a  relative  water  depth 

2iTa/L  relative  wave  length 
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I.   INTRODUCTION 

In  recent  years,  there  has  been  an  accelerated  interest 
in  the  exploration  of  the  ocean  bottom  with  special  atten- 
tion devoted  to  those  depths  encompassed  by  the  Continental 
Shelf.   This  interest  has  been  motivated  in  great  part  by 
the  desire  to  exploit  the  mineral  assets  of  the  ocean  bot- 
tom.  Additionally,  investigation  of  ocean  bottom  pipelines 
as  specifically  applied  to  sewage  outfall  lines,  oil  pipe- 
lines and  large  aqueducts  are  under  consideration. 

As  a  result,  there  has  been  considerable  interest  in 
the  hydrodynamic  interaction  of  gravity  waves  with  large 
submarine  structures  resting  on  the  ocean  bottom.   The 
design  of    large  submarine  offshore  structures  necessarily 
requires  one  to  have  an  understanding  of  the  forces  induced 
by  such  environmental  phenomena  as  surface  waves.   To  this 
point,  the  theoretical  work  has  yet  to  be  a  conclusive  tool 
in  assisting  the  design  engineer  to  predict  all  expected 
conditions  to  which  the  anchored  structure  will  be  subjected 
as  it  lies  on  the  ocean  floor.   Experimental  observations 
of  the  forces  induced  by  ocean  waves  are  difficult  to  ac- 
quire, and  it  is  therefore  necessary  to  simulate  as  pre- 
cisely as  possible  the  conditions  existing  at  the  ocean 
bottom. 

Extensive  experimental  work  has  been  carried  out  since 
Stokes,  who  in  one  of  the  earliest  papers  on  the  subject 
[1],  showed  that,  by  experimentation  with  a  pendulum- type 


11 


system,  the  expression  for  the  force  on  a  cylinder  oscilla- 
ting in  an  infinite  fluid  consists  of  two  terms,  one  in- 
volving acceleration   of  the  cylinder,  and  the  other 
involving  velocity. 

Wave  forces  exerted  on  piles  have  been  intensively  in- 
vestigated during  the  past  two  decades.   It  can  be  shown 
that  when  the  relative  size  of  the  body,  i.e.,  the  cylinder 
diameter,  is  small  compared  to  wave  length,  a  condition 
which  is  generally  satisfied  in  the  case  of  piles,  the 
familiar  Morison  equation  [2]  is  a  valid  approach.   This 
equation  combines  both  a  drag  term  and  an  inertia  term, 
and  the  values  of  the  experimental  coefficients  of  drag 
and  inertia  are  left  to  the  ingenuity  of  the  designer.   At 
best,  these  two  experimental  coefficients  arc  dependent 
upon  the  geometry  of  the  objects  under  consideration  and 
the  amplitude  of  the  fluid  motion,  and  it  would  appear  that 
20  years  of  research  has  shed  little  light  on  understanding 
this  interdependence.   Moreover,  the  proximity  of  the  free 
surface  as  well  as  the  impermeable  bottom  surface  affect 
these  coefficients.   In  view  of  these  difficulties,  the 
Morison  equation  is  not  relied  upon  in  the  present  approach 
to  the  problem,  and  the  results  are  represented  directly  in 
dimensionless  form  rather  than  through  the  drag  and  inertia 
coefficient . 

Keulegan  and  Carpenter  [3]  used  periodic  standing 
waves  in  their  analysis  of  a  cylinder  subjected  to  a  field 
of  oscillatory  motion.   They  were  able  to  separate  the  total 
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resistance  into  drag  and  inertial  components.   They  correla- 
ted values  of  inertia  and  drag  coefficients  under  conditions 
of  varying  current  intensity  and  cylinder  sizes  with  a  peri- 
od parameter,  U  T/D,  a  dimensionless  group  combining  the 
maximum  fluid  velocity  U  ,  the  wave  period  T,  and  the  cylin- 
der diameter.   This  parameter  is  a  form  of  the  relative 
fluid  displacement  since  U  T  is  proportional  to  the  ampli- 
tude of  the  fluid  motion. 

Sarpkaya  and  Garrison  [4]  studied  both  theoretically 
and  experimentally  unidirectional  flow  past  a  circular  cyl- 
inder with  constant  acceleration.   They  also  showed  that 
the  total  resistance  as  well  as  the  drag  and  inertia  com- 
ponents of  force  may  be  represented  as  a  function  of  rela- 
tive fluid  displacement.   Tn  the  present  analysis  of  the 
experimental  data,  this  same  general  kind  of  dependence  on 
fluid  displacement  is  supposed,  and,  accordingly,  the  force 
coefficients  are  correlated  with  the  ratio  of  the  wave 
height  to  cylinder  diameter. 

Johnson  [5]  has  investigated  the  horizontal  forces 
resulting  from  the  interaction  of  gravity  waves  with  a 
circular  cylinder  placed  relatively  near  the  bottom.   Since 
only  long  waves  were  studied,  the  resulting  forces  were 
supposed  to  be  dependent  on  the  water  depth  and  wave  height 
only.   This  study  shed  little  light,  however,  on  the  basic 
problem . 

Herbich  and  Shank,  in  a  recent  study   [6],  have  con- 
ducted wave  channel  testing  of  vertical  and  horizontal  loads 
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on  submerged  three-dimensional  structures  due  to  interaction 
with  gravity  waves.   Several  models  were  studied  using  a 
two-dimensional  wave  tank,  and  they  concluded  that  the  re- 
sults of  these  experimental  studies  were  sufficient  to  be 
applied  to  preliminary  designs.   Unfortunately,  few  actual 
experimental  points  were  presented  and  therefore  the  degree 
of  correlation  is  somewhat  suspect. 

Garrison,  et   al .     [7],  investigated  the  wave/structure 
interaction  problem  as  applicable  to  larger  objects  where 
the  Morison  equation  is  no  longer  valid  and  a  diffraction 
theory  must  be  employed.   In  this  situation,  the  incident 
wave  is  scattered  upon  encountering  the  large  structure, 
and  the  assumption  that  the  object  does  not  affect  the  in- 
cident wave  is  no  longer  valid.   The  diffraction  theory 
accounts  for  the  relative  size  of  the  body  and,  neglecting 
viscous  effects,  assumes  the  existence  of  a  velocity  poten- 
tial.  When  large  structures,  such  as  the  recently  deployed 
Chicago  Bridge  and  Iron  Persian  Gulf  oil  storage  tank,  are 
designed,  the  analysis  must  include  consideration  for  rela- 
tive size  and  free  surface  effect,  and  this  basic  consid- 
eration is  inherent  in  the  diffraction  theory. 

It  appears  that  there  will  continue  to  be  heavy  re- 
liance on  experimental  model  results  in  the  design  process 
until  more  reliable  theoretical  predictions  of  wave  induced 
forces  are  assured.   To  this  time,  large  structures  have 
been  placed  on  the  ocean  floor  only  after  extensive  model 
studies  have  been  carried  out. 
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The  present  study  was  undertaken  to  obtain  an  insight 
into  the  interaction  of  gravity  waves  with  large  submarine 
structures  located  at  or  near  the  ocean  floor.   The  study 
includes  direct  measurement  of  horizontal  and  vertical 
forces  on  the  submerged  cylinder.   The  results  are  pre- 
sented in  dimensionless  form  as  a  function  of  the  relative 
water  depth,  the  relative  wave  length,  and  the  relative 
wave  height. 

The  interest  in  attempting  this  investigation  was 
generated  while  participating  in  numerous  discussions  on 
the  problem  under  the  guidance  of  Dr.  C.  J.  Garrison. 
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II .   THEORETICAL  CONSIDERATIONS 

A  schematic  representation  of  the  problem  under  con- 
sideration is  shown  in  Figure  (1) .   A  rigid  circular  cylin- 
der of  radius  a  submerged  with  center  at  depth  d  in  water 
of  depth  h  is  subjected  to  a  train  of  regular  waves  of  wave 
length  L  and  height  H  progressing  in  the  positive  x-direc- 
tion.   The  primary  interest  lies  in  the  horizontal  and  ver- 
tical components  of  force  induced  by  the  waves  on  the 
circular  cylinder. 

The  complete  theoretical  solution  to  the  present  prob- 
lem is  difficult,  even  assuming  the  flow  to  be  unseparated; 
and,  therefore,  in  order  to  proceed  with  an  experimental 
approach,  it  is  first  necessary  to  apply  dimensional  anal- 
ysis to  the  problem.   It  is  known  a   priori    that  the  maximum 
(horizontal  or  vertical)  force  per  unit  length  of  cylinder 
is  dependent  upon  the  following  variables: 

F  or    F  =    f (h,d,a,L,H,p,g,u)  (1) 

xmax  ymax  v    >>>>>k,  t,  >^j  v    j 

where  h  denotes  the  water  depth,  d  the  depth  of  submergence, 
a  the  cylinder  radius,  L  the  wave  length,  H  the  wave  height, 
p  the  fluid  density,  g  the  acceleration  of  gravity,  and  u 
the  fluid  viscosity.   The  period  of  the  incident  wave  is  not 
required  in  the  above  relationship  because  a  second  rela- 
tionship already  exists  between  the  wave  period  and  the  three 
variables  h,  H,  and  L.   This  is  true  not  only  in  the  case  of 
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cylinder 


\\\\\\\\\\\\\\\\\\W 


channel   floor 


Figure  1.   Definition  Sketch 
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linear  waves,  but  in  the  case  of  large  amplitude  non- linear 
waves  . 

A  dimensional  analysis  of  the  physical  quantities  in- 
dicated in  Equation  (1)  yields  the  following  dimensionless 
variables : 


F         F 

xmax     ymax 
or  — ' 


2?ra 


,  h/a,  d/a,  H/2a, 


Pga' 


Pga' 


DVgha 


(2) 


The  last  term  on  the  right  hand  side  of  Equation  (2)  rep- 
resents the  ratio  of  the  Froude  number  and  the  Reynolds 
number  and,  accordingly,  is  an  indicator  of  the  ratio  of 
viscous  forces  to  gravity  induced  forces.   If  this  number 
is  large,  corresponding  to  a  large  scale  flow  and  fluids 
of  small  viscositv,  it  mav  be  susnected  that  the  gravita- 
tional  forces  may  dominate  the  flow,  and  the  viscous  ef- 
fects may  be  considered  unimportant.   As  a  result,  this 
parameter  may  be  disregarded  and  Equation  (2)  may  be  writ- 
ten as  : 


xmax 


pga' 


or 


ym 


ax 


=  f 


2-rra 


\ 


,  h/a,  d/a,  H/2a 


pga' 


(3) 


A  second  reason  for  the  disregarding  of  this  parameter 
may  also  be  cited.   In  unidirectional  motion  starting  from 
rest  and  moving  with  constant  acceleration,  it  has  been 
shown  by  Sarpkaya  and  Garrison  [4]  that  the  total  force 
exerted  on  a  circular  cylinder  can  be  expressed  as  a  func- 
tion of  the  relative  displacement  of  the  fluid  only.   Al- 
though the  present  problem  involves  oscillatory  as  opposed 
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to  unidirectional  motion,  it  is  nevertheless  suspected  that 
the  same  general  kind  of  dependence  on  relative  displace- 
ment would  result.   In  the  present  problem,  the  amplitude 
of  the  relative  fluid  motion  is  proportional  to  the  wave 
height  to  cylinder  diameter  ratio  and,  therefore,  H/2a 
represents  an  equivalent  parameter.   In  the  experimental 
program,  an  attempt  is  made  to  correlate  the  two  force  co- 
efficients on  the  left  hand  side  of  Equation  (3)  with  the 
dimensionless  ratios  on  the  right  hand  side. 

Although  the  complete  theoretical  solution  to  the  pro- 
posed problem  is  very  complex,  it  is  possible  to  obtain  an 
approximate  solution  if  certain  simplifying  assumptions  are 
applied.   The  velocity  potential  associated  with  the  regular 
progressive  wave  alone  is  given  by: 

g}l     cosh  k(y+h) 

$  =  — r-^rr cos  (kx-at)  (4) 

a    cosh  kh  v     *  K 

and  from  this  expression,  the  velocity  and  acceleration  of 
the  fluid  particles  within  a  wave  may  be  written,  respec- 
tively ,  as  : 

gk^L  cosh  k(y+h) 

u  ="* — -  r  ,  7- —  sin  (kx-at)  (5) 

a     cosh  kh        *•     J 

gk"  sinh  k(y+h) 

v  =  — -  t— vr —  cos  (kx-at)  (6) 

a     cosh  kh        *■     J 

,H 

§K?  cosh  k(y+h)       ri    .>.  (1^ 

u  =  o  cosh  kh    C0S  (kx"at)  (7) 

*  2  sinh  k(v+h)    -   rl    .-,  ,  o-\ 

v  = ut\ — ~  srn  (kx-at)  .  (oj 

a     cosh  kh        ^     J 
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In  wave  force  problems  involving  objects  whose  dimen- 
sions are  small  in  comparison  to  the  wave  length  of  the  in- 
cident wave,  it  is  common  to  assume  that  the  object  does 
not  affect  the  incident  wave.   It  is  then  assumed  that  the 
object  is  subjected  to  a  uniform  oscillatory  flow  having 
acceleration  equal  to  that  induced  by  the  wave  at  the  cen- 
ter of  the  object  if  the  object  were  not  there.   Following 
this  line  of  reasoning,  the  horizontal  force  per  unit  length 
of  cylinder  may  be  written  as  the  product  of  the  virtual 
mass  and  fluid  acceleration  as 


F   =  pira2  (1  +  C  )  Uma,. 


(9) 


if  it  can  be  assumed  that  flow  separation  does  not  occur. 
Substituting  for  u  from  Equation  (7)  and  evaluating  at 
y=d  gives  the  dimensionless  horizontal  force  coefficient 
as 


pga' 


tt(1  +  C  ) 

v    nr 


2i,a)  cosh  T7  (h-d) 
/    cosh  -y  h 


_H 
2a 


(10) 


Since  the  dimensionless  force  coefficient  is  linear  in  the 
wave  height,  this  expression  may  be  written  more  conveniently 
as 


,/,,;0 — ST  =  ff(l+C  ) 

pga2(H/2a)     **    nr 


2fTa 


9-TTO 

cosh  —j —  (h/a  -  d/a) 

,  2?ra  ,  / 
cosh  — f —  h/a 


(ID 


The  parameter  C   occurring  in  Equation  (11)  denotes  the 


m 


added  mass  coefficient  and  depends  on  the  shape  of  the  object 
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in  question  as  well  as  the  proximity  of  the  bottom  and  free 
surface.   Values  for  the  configuration  in  question  which 
corresponds  to  a  circular  cylinder  between  a  rigid  boundary 
and  free  surface  are  not  available  to  this  writer's  know- 
ledge.  The  attempt  in  this  thesis,  therefore,  is  to  deter- 
mine by  trial  and  error  the  value  of  C  which  gives  the 
best  fit  to  the  experimental  data. 
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III.   EXPERIMENTAL  INVESTIGATION 

A.    APPARATUS 

The  wave  channel  as  depicted  in  Figure  (2)  was  de- 
signed with  a  rectangular  cross  section  two  feet  deep  with 
a  width  of  15  inches.   The  channel  was  fabricated  in  sec- 
tions.  Each  section  was  composed  of  two  vertical  walls 
and  one  horizontal  floor  member.   The  vertical  walls  were 
stiffened  by  use  of  three  4  foot  long  2  x  4's  attached  to 
3/4  inch  x  8  foot  x  3  foot  sheets  of  plywood  at  the  ends 
and  middle  of  its  longest  dimension.   The  channel  floor 
member  was  composed  of  two  15  inch  x  8  foot  long  plywood 
sheets  separated  by  a  1  inch  x  4  inch  strip  around  the 
edge.   The  bottom  and  side  members  were  mated  together 
with  long  3/8  inch  threaded  rods  and  formed  the  basic 
channel  cross  section.   Six  of  these  8  foot  sections  were 
then  bolted  together  to  form  the  basic  channel.   Square- 
ness of  the  channel  cross  section  was  maintained  by  use 
of  spacers  placed  between  the  stiffening  members  near  the 
floor.   Every  butt  joint  was  given  a  heavy  application  of 
silicone  sealant  to  prevent  leakage. 

At  a  location  30  feet  from  the  wave  generator  end  of 
the  channel,  a  3  foot  wide  plexiglass  observation  window 
extending  the  full  height  of  the  channel  was  installed  in 
both  walls  of  the  channel  as  shown  in  Figure  (3) .   The 
test  cylinder  was  then  mounted  in  the  center  of  this  win- 
dow. 
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Figure  3.   Wave  Channel  Observation  Window 
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Dissipation  of  the  generated  waves  was  provided  by  the 
beach  structure  located  in  the  last  section  of  the  channel. 
The  beach  as  shown  in  Figure  (4)  consisted  of  three  layers 
of  perforated  stainless  steel  sheeting.   Layer  separation 
and  structural  framing  were  provided  by  2  x  4  runners  with 
their  largest  cross  section  dimension  perpendicular  to  the 
sheeting.   The  slope  of  the  channel  beach  was  maintained 
at  6:1. 

The  wave  generating  component  of  the  channel  as  shown 
in  the  photograph  of  Figure  (5)  had  two  basic  elements,  the 
piston-type  plate  and  the  driving  motor.   The  plate  was  an 
aluminum  sheet  1/8  inch  thick  with  adjustable  edges  on  the 
vertical  sides  giving  it  the  capability  of  adjusting  to 

LOiuOii.i     lO      cue     dxCii^o     w  x      tlic     s-.li  ciilitw  x  .  j.  iic     pj.uco      v\ai     iuuuhucu 

on  a  framework  attached  to  4  slider  bearings,  and  the 
slider  bearings  rode  on  two  1.0  inch  brass  rods  attached 
to  the  channel  walls.   A  driving  rod  was  attached  to  the 
plate  structural  members  and  connected  to  a  flywheel  at- 
tached to  the  output  shaft  of  a  Vari-Drive  motor.   The 
drive  rod  was  attached  to  the  flywheel  by  means  of  a  bolt 
and  slot  in  such  a  manner  that  the  eccentricity  could  be 
varied  between  zero  and  6  inches.   In  this  way,  waves  of 
various  amplitudes  could  be  generated.   The  variable  speed 
transmission  was  driven  by  a  2  hp  inductance  motor  and  pro- 
vided a  speed  range  extending  from  22  to  180  rpm  giving 
wave  lengths  ranging  from  18  feet  to  1.5  feet  at  a  water 
depth  of  18  inches. 
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A  plexiglass  circular  cylinder  of  6  inch  diameter  with 
an  end  clearance  of  approximately  1/16  inch  was  used  in  the 
test  reported  herein.   As  depicted  in  Figure  (6) ,  the  cylin 
der  was  hung  by  two  fine  wires  at  the  appropriate  distance 
off  the  wave  channel  floor.   Weights  were  added  inside  the 
cylinder  to  pretension  the  wires  and  to  prevent  slack  in 
the  wires  during  testing. 

The  upper  ends  of  the  wires  were  attached  to  small 
aluminum  cantilever  beams  equipped  with  strain  gages  on 
both  the  compression  and  tension  side  of  the  beams.   These 
vertical  force  strain  gages  were  then  connected  into  a 
Wheatstone  bridge  and  read-out  on  one  channel  of  a  two- 
channel  Sanborn  Model  296  amplifier/recorder.   The  cylinder 
was  restrained  in  the  horizontal  direction  by  similar  fine 
wires  connected  into  a  harness  and  around  a  4  inch  diameter 
ball  bearing  pulley  mounted  on  the  wave  channel  floor  as 
shown  in  Figure  (7).   These  wires  were  then  attached  to  a 
pair  of  cantilever  beams  equipped  with  strain  gages,  and 
the  strain  gages  were  connected  into  a  Wheatstone  bridge 
and  read-out  on  the  second  channel  of  the  two- channel  am- 
plifier/recorder. 

The  four  beams,  two  for  the  vertical  force  and  two  for 
the  horizontal  force,  were  attached  to  a  plywood  mounting 
board  positioned  over  the  channel  as  shown  in  Figure  (8). 
The  beams  were  all  exactly  alike,  and  the  section  of  the 
beams  where  the  strain  gages  were  mounted  was  reduced  to  a 
cross  section  of  1/4  x  1/4  inches  from  a  nominal  beam  cross 
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Figure  7.   Horizontal  Wire  Harness  and  Pulley  Assembly 
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Figure  8.   Beam  Mounting  Board. 
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section  of  1/2  x  3/8  inches.   The  machining  provided  a 
smooth  transition  into  the  area  of  stress  concentration. 
SR-4  FAE-25-12S13  strain  gages  were  applied  to  the  top 
and  bottom  surfaces  of  the  beams.   Each  beam  was  9  inches 
long  and  had  mounting  holes  for  two  moment  arm  positions, 
one  at  6  inches  and  one  at  3  inches  from  the  strain  gage 
location.   This  variable  load  attachment  point  gave  each 
beam  a  sensitive  force  measurement  capability  under  a 
large  range  of  values.   The  beams  were  designed  to  with- 
stand a  six-pound  load  at  the  largest  moment  arm. 

The  wave  height  gage,  as  shown  in  Figure  (9) ,  con- 
sisted of  a  small  diameter  aluminum  rod  acting  as  the 
basic  frame  to  which  a  plexiglass  mounting  table  was  at- 
tached.  The  elements  of  the  Wheatstone  bridae  were 
mounted  on  this  table,  and  two  30-gage  copper  wires  spaced 
1/2  inch  apart  were  extended  from  the  table  and  attached 
to  an  insulator  mounted  at  the  end  of  the  aluminum  rod. 
An  18-inch  separation  between  the  levels  of  the  table  and 
insulator  provides  a  large  measurement  range  capability. 
The  other  end  of  the  rod  was  mated  to  a  traverse  block 
stand  which  allowed  the  gage  to  travel  ±  6  inches,  and 
provided  a  convenient,  yet  accurate,  means  for  calibration 
by  adjusting  the  submergence  of  the  wires.   The  complete 
unit  was  positioned  over  the  channel  and  centered  at  ap- 
proximately 5  feet  in  front  of  the  cylinder. 

The  wave  height  probe  was  of  the  parallel  wire  re- 
sistance type.   This  type  of  gage  operates  on  the  principle 
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Figure  9.   Wave  Height  Gage. 
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that  the  conductance  between  the  two  parallel  submerged 
wires  varies  approximately  proportionately  to  the  length 
of  wire  submerged.   The  wave  gage  was  connected  in  parallel 
to  one  leg  of  a  Wheatstone  bridge  as  depicted  in  Figure 
(10) .   To  compensate  for  the  inability  of  the  carrier  am- 
plifier to  offset  the  resistance  of  the  immersed  parallel 
wires,  an  external  variable  10K  ohm  potentiometer  was  con- 
nected across  an  opposite  leg  of  the  bridge.   This  resistor 
was  used  for  coarse  balancing  of  the  bridge.   The  Wheatstone 
bridge  was  connected  via  signal  and  excitation  leads  to  a 
portable  carrier  preamplifier  and  read  out  on  a  one-channel 
H/P  Sanborn  301  amplifier/recorder. 

During  preliminary  operation  while  attempting  to  clear 
the  problems  from  the  system,  it  was  found  that  installation 
of  an  impermeable  wall  between  cylinder  and  channel  floor 
was  a  necessity.   The  experiment  was  designed  to  simulate 
the  placement  of  the  circular  cylinder  on  the  bottom,  but 
it  was  necessary  to  maintain  a  small  but  finite  gap  between 
cylinder   and  the  bed  of  the  channel  to  insure  that  the 
pair  of  wires  connected  to  the  cantilever  beams  provided 
the  sole  vertical  support.   However,  it  was  found  that  the 
vertical  forces  as  well  as,  to  a  lesser  extent,  the  hori- 
zontal forces  were  strongly  dependent  on  the  gap  width. 
Flow  visualization  studies  showed  that  a  high  velocity 
jet  occurred  in  the  small  gap  under  the  cylinder.   Thus, 
in  order  to  eliminate  the  flow  under  the  cylinder  and  yet 
not  influence  the  force  measurement  system,  a  flexible 
plastic  sheet  was  installed  across  the  gap  between  the 
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Figure  10.   Wave  Height  Gage  Bridge  Schematic 
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cylinder  and  floor.   The  sheet  was  joined  to  the  cylinder 
by  pressing  the  plastic  into  a  machined  slot  in  the  cylin- 
der wall,  and  then  anchoring  it  with  a  length  of  "0"-ring. 
The  other  end  of  the  plastic  was  laid  under  a  metal  sheet 
12  inches  x  14  1/2  inches  and  1/16  inch  thick  and  placed 
at  the  channel  bottom  as  shown  in  Figure  (11)  .   The  gap 
between  the  cylinder  and  floor  was  adjusted  to  approxi- 
mately 1/8  inch.   This  gap  was  found  to  be  large  enough 
to  insure  that  the  cylinder  never  touched  the  bottom  when 
subjected  to  the  largest  wave. 

In  order  to  demonstrate  the  rather  large  influence  of 
this  small  gap  under  the  cylinder,  Figure  (12)  was  prepared 
showing  both  the  maximum  horizontal  and  vertical  force  as 
a  function  of  wave  amplitude.   It  is  apparent  from  this 
figure  that  removal  of  the  plastic  barrier  had  the  effect 
of  reducing  the  maximum  horizontal  force  a  slight  amount. 
With  the  barrier  in  place,  the  upward  force  was  generally 
quite  large  while  the  downward  force  was  rather  small. 
However,  upon  removal  of  the  barrier,  a  rather  high  veloc- 
ity flow  through  the  gap  was  observed  by  injection  of  dye 
near  the  wave  channel  floor,  and  this  local  region  of  high 
velocity  apparently  resulted  in  low  pressures  on  the  under- 
side of  the  cylinder  causing  the  downward  force  to  be  in- 
creased considerably.   The  rather  extreme  effect  of  the 
gap  is  clearly  demonstrated  by  comparison  of  the  various 
curves  indicated  on  Figure  (12) . 
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Inasmuch  as  the  intent  of  the  present  study  was  to 
determine  wave  forces  on  a  cylinder  lying  on  the  bottom 
wherein  the  gap  would  not  be  present,  the  plastic  barrier 
was  maintained  in  place  throughout  the  test  program  for 
all  cases  where  the  cylinder  was  placed  near  the  bottom. 

B.    TEST  PROCEDURE 

The  first  step  in  the  test  procedure  was  to  fill  the 
wave  channel  with  water  to  obtain  the  desired  water  depth 
to  cylinder  diameter  ratio,  h/a.   Experiments  were  con- 
ducted at  the  three  depths,  12,  15,  and  18  inches.   After 
the  desired  water  level  was  attained,  a  careful  inspection 
was  carried  out  to  insure  that  all  equipments  were  properly 
prepared  for  operation. 

The  vertical  wire  position  screws  were  adjusted  to 
obtain  the  desired  elevation  of  the  cylinder  off  the  chan- 
nel bottom.   The  horizontal  wires  were  pretensioned  to 
prevent  slacking  and  the  attendant  errors  in  force  measure- 
ments . 

A  typical  run  may  be  described  as  follows.   With  water 
adjusted  to  the  proper  depth,  a  given  wave  length,  or  equiv- 
alently,  wave  period  was  established  by  setting  the  speed 
of  the  Vari-Drive  unit.   Subsequent  to  a  sufficient  warm-up 
period,  the  amplifier/recorders  shown  in  Figure  (13)  were 
balanced  and  readied  for  calibration.   Calibration  of  the 
three  recording  channels,  horizontal  force,  vertical  force, 
and  wave  height  was  carried  out  at  the  attenuation  expected. 
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to  be  in  use  during  the  run.   After  calibration,  the  wave 
height  was  set  by  adjusting  the  eccentricity  of  the  wave 
generator  to  the  appropriate  amplitude.   The  wave  generator 
was  then  turned  on,  and  after  a  steady  state  wave  was  es- 
tablished in  the  channel,  the  recorders  were  turned  on  and 
data  was  taken.   The  wave  generator  was  then  switched  off, 
and  the  next  wave  height  was  set  on  the  flywheel  eccentric 
while  maintaining  the  same  speed  setting.   Wave  heights  at 
the  fixed  wave  length  were  varied  within  the  capability  of 
the  system  to  produce  representative  readings  usually  ranging 
from  approximately  1/4  inch  to  5  inches.   At  the  conclusion 
of  each  series  of  runs  at  a  given  wave  length,  a  check  cali- 
bration on  the  three  channels  was  quickly  carried  out.   Then, 
a  new  wave  length  was  selected  and  set  on  the  Vari-Drive, 
and  the  cycle  of  balancing,  calibration,  wave  height  selec- 
tion, and  wave  production  was  repeated. 

Calibration  was  performed  on  the  load  cells  as  follows. 
The  recording  chart  was  set  at  1   /sec.  speed,  and  loads  in 
one-pound  steps  were  placed  on  the  saddles  located  beneath 
and  hung  from  the  horizontal  sensing  wires  as  shown  in  Fig- 
ure (14) .   Calibration  loads  were  placed  initially  on  the 
after  beam,  then  on  the  forward  beam.   The  vertical  load 
cell  was  calibrated  by  use  of  a  platform  of  known  weight 
centered  on  the  cylinder  as  shown  in  Figure  (15).   Weights 
were  then  placed  on  this  platform. 

Wave  height  gage  calibration  was  performed  by  adjusting 
the  vertical  position  of  the  gage  through  the  use  of 
bearing  and  screw.   The  procedure  was  carried  out  as  follows. 
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The  gage  was  raised  1/2  inch,  simulating  a  decrease  in  water 
level,  then  lowered  1/2  inch  below  the  null,  simulating  an 
increase  in  water  level.   Then,  the  gage  was  raised  +1  inch 
and  lowered  -1  inch,  etc.   This  procedure  closely  simulated 
the  actual  passage  of  a  wave  profile  with  crests  and  troughs 
The  generated  wave  crest  passing  through  the  test  location 
is  shown  in  Figure  (11)  . 
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IV.   EXPERIMENTAL  RESULTS 

As  indicated  in  Equation  (3)  in  the  THEORETICAL  CON- 
SIDERATIONS, the  wave  force  coefficients  can  be  represented 
as  functions  of  relative  water  depth,  relative  depth  of  sub- 
mergence, relative  wave  length,  and  relative  wave  height  pa- 
rameters.  Accordingly,  a  test  program  was  devised  to  measure 
the  wave  force  coefficients  for  practical  ranges  of  the  in- 
dependent parameters.   Four  series  of  tests  were  conducted 
at  different  water  depths  and  depths  of  submergence.   Three 
of  these  corresponding  to  depth  to  cylinder  radius  ratios 
of  4.0,  5.0,  and  6.0,  were  carried  out  with  the  cylinder 
near  the  bottom,  i.e.,  at  d/a=  3.0,  4.0,  and  5.0,  respec- 
tively.  A  fourth  series  of  tests  was  conducted  with  the 
cylinder  raised  one  radii  off  the  bottom  corresponding  to 
d/a=  4.0  at  a  water  depth  ratio  of  6.0. 

A.    TYPICAL  WAVE  FORCE  TRACES 

Inasmuch  as  it  is  difficult  to  completely  characterize 
the  force  variations  with  a  single  number,  typical  strip 
chart  read-outs  for  the  various  series  of  tests  are  pre- 
sented in  Figures  (16-28).   Figures  (16-19)  show  the  hori- 
zontal and  vertical  force  variations  with  the  cylinder  on 
the  bottom  corresponding  to  a  relative  depth  of  submergence 
of  d/a  =  5.0  and  a  relative  water  depth  ratio  of  h/a=  6.0. 
The  figures  are  arranged  according  to  increasing  values  of 
the  wave  length  parameter  ranging  from  2Tra/L=  0.09  to  0.97. 
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On  each  figure,  data  corresponding  to  four  values  of  rela- 
tive wave  height  are  presented  in  ascending  order. 

Figure  (16)  representing  the  largest  wave  length  and 
having  a  depth  to  wave  length  ratio  of  h/L=  0.08  corre- 
sponds to  waves  which  are  somewhat  outside  the  sinusoidal 
range.   These  long  waves  are  characterized  by  long,  smooth 
troughs  and  sharp  crests.   Accordingly,  the  horizontal 
force  traces  are  characterized  by  positive  and  negative 
peaks  which  result  from  the  maximum  acceleration  occurring 
on  either  side  of  the  crest.   These  two  peaks  are  followed 
by  a  relatively  flat  portion  of  the  force  trace  corres- 
ponding to  the  long,  smooth  trough  portion  of  the  wave.   As 
can  be  seen  from  the  figure,  this  kind  of  variation  appears 
to  predominate  even  as  the  wave  height  is  increased  for 
this  long  wave  case. 

The  vertical  force,  however,  shows  a  much  more  incon- 
sistent variation.   If  the  system  natural  frequency  is  dis- 
regarded, definite  upward  peaks  are  observed  at  points  on 
the  trace  corresponding  to  the  wave  crests.   These  peaks 
are  apparently  related  to  the  reduced  pressure  over  the 
upper  portion  of  the  cylinder ' which  is  associated  with  the 
velocity  squared  term  in  the  Bernoulli  equation  since  the 
maximum  velocity  occurs  at  the  crests.   In  sinusoidal  waves, 
an  equal  and  opposite  velocity  occurs  at  the  troughs  which 
would  naturally  give  rise  to  a  second  upward  peak  in  the 
force  trace.   However,  in  shallow  water  waves,  the  veloc- 
ity associated  with  a  trough  is  much  reduced   and,  there- 
fore peaks  corresponding  to  the  crests  only  are  observed. 
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Figure    16.      Wave    Force   Trace,    h/a=    6.0,    d/a=    5.0,    2tTa/L=    0.09 
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Figure  (17)  corresponds  to  a  shorter  wave  length  having 
a  relative  wave  length  value  of  2iTa/L=  0.16.   It  is  apparent 
from  the  figure  that  both  the  horizontal  and  vertical  forces 
vary  sinusoidally  for  the  smallest  amplitude  wave  correspond- 
ing to  H/2a=  0.5.   As  the  wave  height  is  increased,  it  is  ap- 
parent that  the  horizontal  force  continues  to  vary  sinusoid- 
ally  for  the  full  range  of  amplitudes.   The  vertical  force, 
however,  shows  a  sinusoidal  variation  for  the  smallest  am- 
plitude waves  and  becomes  increasingly  less  sinusoidal  as 
the  wave  height  is  increased.   Moreover,  two  separate  peaks 
appear  to  develop  on  the  vertical  force  trace  as  wave  height 
increases,  one  associated  with  the  maximum  velocity  occur- 
ring at  the  crest  and  the  other  associated  with  the  trough. 
The  maximum  velocity  at  both  the  crest  and  the  trough  re- 
sults in  reduced  pressures  on  the  upper  portion  of  the  cylin- 
der and,  consequently,  corresponds  to  upward  vertical  forces 

Figures  (18)  and  (19)  correspond  to  still  smaller  wave 
lengths  of  27ra/L=  0.24  and  0.77,  respectively.   It  may  be 
noted  that  again  the  horizontal  force  shows  a  sinusoidal 
variation  throughout  the  range  of  wave  amplitudes  while  the 
vertical  force  shows  the  sinusoidal  variation  at  small  am- 
plitudes becoming  non-sinusoidal  as  the  amplitude  is  in- 
creased.  Figure  (19),  however,  which  corresponds  to  the 
shortest  wave  lengths  shows  sinusoidal  variations  through- 
out the  range  of  wave  amplitudes  for  the  horizontal  as  well 
as  vertical  force. 
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igure  17.   Wave  Force  Trace,  h/a=  6.0,  d/a=  5.0,  2-a/L=  0.16 
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Figure    18.      Wave    Force   Trace,    h/a=    6.0,    d/a=    5.0,    2ira/L= 
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Figure    19.      Wave    Force   Trace,    h/a=    6.0,    d/a=    5.0    2ira/L  =    0.7  7 
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Horizontal  and  vertical  force  coefficients  correspond- 
ing to  relative  water  depth  of  h/a=  5.0  with  the  cylinder 
on  the  bottom,  i.e.,  d/a=  4.0  are  shown  in  Figures  (20-22). 
These  three  figures  represent  three  different  values  of 
relative  wave  length  extending  from  27Ta/L=  0.10  to  0.67. 
The  general  trends  tend  to  be  similar  to  those  corresponding 
to  the  larger  depth  of  h/a  as  shown  in  Figures  (16-19)  and 
previously  discussed.   Data  corresponding  to  the  smallest 
value  of  2iTa/L  as  presented  in  Figure  (20)  shows  trends 
typical  of  nonlinear  shallow  water  v/aves.   The  horizontal 
force  shows  the  null  period  associated  with  the  long,  smooth 
troughs,  while  the  vertical  force  is  characterized  by  the 
upward  force  peak  corresponding  to  the  wave  crests.   As  the 

1- T  D  "\  r  *"*        1   ^  "P  T  "*~  V*        1    C        Ar\r*Trr\ryc~r\A  4-  V>  r\       V>nr  i    nn-nt  n  1         -Pa  vrn        n  rv  1  t  Tt        rinc- 

plays  sinusoidal  variation  as  indicated  in  Figures  (21)  and 
(22).   For  the  shorter  wave  lengths,  the  vertical  force  shows 
a  sinusoidal  variation  at  small  amplitudes  with  superimposed 
upward  peaks  occurring  near  the  maximum  velocity  location  at 
the  crests  and  troughs.   Figure  (22)  corresponding  to  the 
shortest  wave  length  shows  a  sinusoidal  variation  for  both 
force  components  throughout  the  range  of  wave  heights. 

Typical  force  variation  records  for  the  series  of  tests 
conducted  at  h/a=  4.0  are  presented  in  Figures  (23-25). 
Three  representative  values  of  relative  wave  length,  2iTa/L  = 
0.11,  0.21,   and  .45  are  presented.   Figure  (23)  shows  varia- 
tions similar  to  those  of  equal  values  of  2-rra/L  at  the  two 
larger  depths  already  discussed.   Figures  (24)  and  (25) 
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Figure    20.      Wave    Force   Trace,    h/a=    5.0,    d/a=    4.0,    2Tra/L=    0 .  10 
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Figure    21.      Wave    Force   Trace,    h/a=    5.0,    d/a=    4.0,    2iTa/L=    0.19 
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Figure    22.      Wave    Force   Trace,    h/a=    5.0,    d/a=    4.0,    2TTa/J 
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Figure    23.      Wave    Force   Trace,    h/a=    4.0.    d/a=    3.0,    2ira/L=    0.11 
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also  similar  to  previous  runs,  the  primary  difference  being 
the  increased  magnitude  occurring  at  the  smaller  depths. 

The  final  series  of  runs  was  carried  out  at  a  relative 
water  depth  ratio  of  h/a=  6.0  with  the  cylinder  suspended 
one  radii  off  the  bottom,  i.e.,  d/a=  4.0.   Typical  force 
traces  at  three  different  values  of  relative  wave  length, 
2iTa/L=  0.10,  0.16,  and  0.6,  are  presented  in  Figures  (26), 
(27),  and  (28),  respectively.   The  horizontal  force  trace 
as  shown  on  the  figures  reflects  the  same  general  trends  as 
the  previously  described  series.   The  vertical  component 
of  force  also  shows  trends  similar  to  those  previously  dis- 
cussed with  the  cylinder  on  the  bottom,  but  the  peaks  as- 
sociated with  maximum  velocity  in  the  wave  occurring  at  the 
crests  and  troughs  are  seen  to  be  in  the  opposite  direction. 
That  is,  the  velocity  of  the  fluid  through  the  one  radii  gap 
between  the  cylinder  and  the  floor  is  apparently  consider- 
ably larger  than  the  velocity  through  the  space  above  the 
cylinder  and  beneath  the  free  surface  and,  therefore,  a 
sucking-down  action  as  opposed  to  a  lifting  action  results. 
As  in  previous  runs,  it  may  be  noted  that  Figure  (28),  cor- 
responding to  the  largest  value  of  2iTa/L,  shows  a  sinusoidal 
variation  for  both  force  components  throughout  the  range  of 
wave  heights . 

B.    WAVE  FORCE  COEFFICIENTS 

In  order  to  represent  the  results  of  the  experiments  in 
a  simpler,  albeit  less  descriptive  form,  the  maximum  values 
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Figure    24.      Wave    Force   Trace,    h/a=    4.0,    d/a=    3.0,    2Tra/L=    0.21 
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Figure  25.   Wave  Force  Trace,  h/a=  4.0,  d/a=  3.0,  2tt. 


59 


+ 


*\y 


Figure  26.   Wave  Force  Trace,  h/a=  6.0,  d/a=  4.0,  2      O.K 
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Figure    27.      Wave    Force   Trace,    h/a=    6.0,    d/a=    4.0,    27Ta/L  =    0. 
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Figure    28.      Wave    Force   Trace,    h/a=    6.0,    d/a=    4.0,    2ira/L=    0.32 
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of  the  dimensionless  force  coefficients  obtained  from  the 
force  traces  were  plotted  as  a  function  of  relative  wave 
height,  H/2a.   Results  corresponding  to  the  relative  water 
depth  of  h/a^  6.0  with  the  cylinder  on  the  bottom,  i.e., 
d/a=  5.0,  are  presented  in  Figures  (29-32).   Results  cor- 
responding to  the  other  two  relative  water  depths  of  h/a= 
5.0  and  4.0  also  with  the  cylinder  on  the  bottom  are  pre- 
sented in  Figures  (33-36)  and  (37-40),  respectively.   These 
results  correspond  to  values  of  2iTa/L  ranging  from  approxi- 
mately 27Ta/L=  0.08  to  0.97. 

As  was  apparent  from  the  force  traces,  the  horizontal 
force  showed  sinusoidal  variations  with  equal  magnitude  in 
both  directions  and,  therefore,  only  one  curve  was  required 
to  represent  the  amplitude  of  this  variation  with  relative 
wave  height.   The  dimensionless  force  coefficients  are  de- 
fined as  f=F/pga2 ,  and  three  curves  are  noted  on  each  fig- 
ure, one  for  the  horizontal  and  one  each  for  the  upward  and 
downward  vertical  force.   As  is  evident  from  the  figures, 
the  horizontal  force  coefficient  shows  a  linear  variation 
with  the  wave  height  in  all  cases. 

However,  as  the  vertical  force  showed,  in  general,  larger 
upward  than  downward  magnitudes,  the  two  curves  labeled  f 
and  f   are  presented  showing  the  maximum  upward  and  downward 
values  of  the  vertical  force  coefficient,  respectively. 
Throughout  the  range  of  the  data,  the  vertical  force  coeffi- 
cient for  the  cylinder  located  on  the  bottom  showed  a  quite 
consistent  variation,  the  upward  force  being  always  larger 
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RELATIVE  WAVE  HEIGHT,  H/2a 
Figure  29.   Wave  Force  Coefficients  for  h/a=  6.0,  d/a=  5.0 
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RELATIVE   WAVE   HEIGHT,    H/2a 
Figure    30.      Wave    Force    Coefficients    for   h/a=    6.0,    .: 


65 


•1.0r 


.8 


.6 


4 


m 


CO 

H 

w 

(—1 
CJ 
(— I 

Pi 

w 
o 
u 

w 
u 
cd 
o 
til 


0 
1.0 


.8 


.6 


27ra/L  =  0.24 


2Tra/L  =  0.49 


f 


x 


2 


0 

RELATIVE   WAVE   HEIGHT,    11/ 2 a 
Figure    31.      Wave    Force    Coefficients    for  h/a=    6.0,    d/a=    5.0 
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Figure  32.   Wave  Force  Coefficients  for  h/a=  6.0,  d/a 
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than  the  downward  force  and  having  a  generally  concave  up- 
ward trend.   The  downward  force  was  always  quite  small  show- 
ing a  nearly  linearly  increasing  variation  followed  by  a 
breaking  off  at  higher  wave  amplitudes. 

The  variations  of  the  force  coefficients  as  observed 
in  Figures  (29-40)  can  be  justified  if  one  considers  the 
components  which  make  up  the  dynamic  pressure  acting  on  the 
surface  of  the  cylinder.   Two  components  of  pressure  are 
generally  recognized,  the  one  component  being  associated 
with  the  unsteady  term  in  the  Bernoulli  equation,  and  the 
second  being  associated  with  the  velocity  squared  term  oc- 
curring therein.   The  unsteady  term  involves  the  velocity 
potential  to  the  first  power  and,  therefore,  is  linear  in 
H/2a.   The  second  term  involves  Lhe  spatial  derivative  of 
the  velocity  potential  squared  and  is,  therefore,  propor- 
tional to  (H/2a)2.   Thus,  linear  variation  of  the  horizon- 
tal force  with  H/2a  apparently  results  from  contributions 
from  the  unsteady  component  of  pressure  only.   The  velocity 
squared  term  which  would  contribute  a  nonlinear  component 
of  force  can  contribute  nothing  to  the  horizontal  force  if 
the  flow  pattern  is  symmetrical  about  the  y-axis.   This, 
of  course,  is  expected  provided  separation  does  not  occur 
and  accounts  for  the  linear  variation  in  the  horizontal 
force  with  the  parameter  H/2a. 

The  vertical  force  variation  with  H/2a,  however,  shows 
markedly  different  trends.   Both  the  upward  and  downward 
components  seem  to  vary  linearly  with  relative  wave  hei 
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RELATIVE  WAVE  HEIGHT,  H/2a 
Figure  53.   Wave  Force  Coefficients  for  h/a=  5.0,  d/a=  4.0. 
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RELATIVE  WAVE  HEIGHT,  H/2a 
Figure  34.   Wave  Force  Coefficients  for  h/a=  5.0,  d/a=  4.0 
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RELATIVE  WAVE  HEIGHT,  H/2a 
Figure  37.   Wave  Force  Coefficients  for  h/a=  4.0,  d/a=  3.0. 
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RELATIVE   WAVE   HEIGHT,    H/2a 
Figure    38.      Wave    Force    Coefficients    for   h/a=    4.0,    d 
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Figure  39.   Wave  Force  Coefficients  for  h/a=  4.0,  d/a=  3. 
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RELATIVE  WAVE  HEIGHT,  H/2a 
Figure  40.   Wave  Force  Coefficients  for  h/a=  4.0,  d/a=  3.0. 
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for  small  amplitude  waves,  and  the  upward  force  tends  to 
bend  upwards  while  the  downward  force  tends  to  decrease  at 
larger  values  of  H/2a.   These  variations,  however,  are  con- 
sistent with  the  previous  discussion  regarding  the  pressure 
on  the  surface  of  the  cylinder.   In  the  case  of  the  vertical 
force,  the  flow  pattern  is  not  symmetrical  about  a  horizon- 
tal axis  passing  through  the  center  of  the  cylinder,  and, 
therefore,  the  contribution  from  the  velocity  squared  term 
in  the  Bernoulli  equation  does  not  cancel.   As  a  crest  or 
trough  passes  the  cylinder,  the  wave  induced  velocity  is  a 
maximum  giving  rise  to  a  component  of  pressure  (or  force) 
proportional  to  (H/2a)2.   The  vertical  force  coefficient 
•variation  with  H/2a,  therefore,  is  composed  of  two  terms, 
one  proportional  to  H/2a  and  one  proportional  to  (H/2a)2. 
Thus,  its  somewhat  parabolic  shape  is  not  unexpected. 

The  maximum  values  of  the  horizontal  and  vertical  force 
coefficients  for  the  final  series  of  runs  conducted  at  a 
water  depth  ratio  of  h/a=  6.0  with  the  cylinder  suspended 
off  the  bottom,  i.e.,  d/a=  4.0,  are  presented  in  Figures 
(41)  through  (44)..  As  in  the  previously  described  tests, 
the  horizontal  force  coefficients  showed  a  linear  variation 
with  wave  amplitude  whereas  the  vertical  force  coefficients 
showed  an  initial  linear  portion  followed  by  a  nonlinear  por 
tion.   However,  in  this  series  of  runs,  the  high  velocity 
in  the  gap  beneath  the  cylinder  caused  a  net  nonlinear  con- 
tribution to  the  force  in  the  downward  direction  making  the 
magnitude  of  f"  larger  than  f  ,  a  trend  opposite  to  those 
associated  with  the  cylinder  near  the  bottom. 
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Figure  41.   Wave  Force  Coefficients  for  h/a=  6.0,  d,     .0. 
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Figure  42.   Wave  Force  Coefficients  for  h/a=  6.0,  d/a=  4.0 
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Figure  44.   Wave  Force  Coefficients  for  h/a=  6.0,  d/a=  4.0. 
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C.    HORIZONTAL  FORCE  COEFFICIENT 

The  maximum  values  of  the  horizontal  force  coefficient 
f  ,  as  shown  in  Figures (29 -40) ,  are  seen  to  vary  linearly 
with  the  relative  wave  height,  H/2a.   In  view  of  this  lin- 
ear behavior,  it  is  possible  to  represent  these  forces  in 
a  much  more  condensed  form.   That  is,  since  the  force  co- 
efficients are  linear  with  H/2a,  the  parameter, 

Fx/[pga2(H/2a)], 

which  actually  represents  the  slopes  of  these  curves  can 
be  presented  in  a  single  figure  as  a  function  of  the  rela- 
tive wave  length,  2-rra/L.   Such  a  plot  is  presented  in  Fig- 
ure (45)  wherein  this  force  coefficient  is  represented  as 
a  function  of  2na/L  for  the  four  series  of  tests  carried 
out  at  different  water  depths  and  depths  of  submergence. 

The  upper  three  curves  shown  in  Figure  (45)  represent 
three  different  water  depths  wherein  the  cylinder  was  lo- 
cated on  the  bottom.   The  effect  of  water  depth  is  clearly 
indicated,  as  the  curves  corresponding  to  decreasing  depths 
lie  above  the  one  at  a  maximum  water  depth  of  h/a=  6.0. 
The  general  trends  of  these  three  curves  are  all  similar, 
maximum  values  of  force  coefficients  occurring  near  2TTa/L  = 
0.3  with  decreasing  force  at  shorter  wave  lengths.   At 
longer  wave  lengths  than  those  corresponding  to  the  maximum 
value  for  the  force  coefficient,  the  curves  tend  to  become 
flat  while  the  linear  theory  would  indicate  a  decrease 
similar  to  that  associated  with  the  curve  correspondin; 
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EXPERIMENTAL       APPROXIMATE  ANALYSIS 

EQUATION  (11) 
h/a     d/a      Cm 


Cx) 


eg 

03 
GO 
Q. 


Uh 


W 
i— i 

u 

I— I 

Uh 
Uh 

w 
o 
c_> 

w 
u 

Pi 
o 


.04  .06  .08  0.1 


Relative  Wave  Length,  2-na/L 
Figure  45.   Horizontal  Force  Coefficients 


0.6  0.8  1.0 


83 


h/a=  6.0  and  d/a=  4.0.   It  is  possible  that  this  behavior 
is  caused  by  the  fact  that  the  wave  is  becoming  a  nonlinear 
shallow  water  wave. 

The  horizontal  force  coefficients  corresponding  to  the 
fourth  series  of  tests  carried  out  at  a  water  depth  ratio 
of  h/a=  6.0  and  with  the  cylinder  suspended  at  one  radii 
off  the  bottom  is  also  shown  in  Figure  (45)  as  the  lowest 
curve.   The  general  trends  are  seen  to  be  somewhat  similar 
to  the  cases  where  the  cylinder  was  located  on  the  bottom 
except  that  the  curve  tends  to  be  somewhat  flatter. 

In  an  attempt  to  compare  the  experimental  results  with 
the  approximate  theory,  Equation  (11)  was  programmed  on  a 
digital  computer,  and  values  for  the  horizontal  force  co- 
efficients were  generated  as  a  function  cf  2ira/L  assuming 

various  values  for  the  added  mass  coefficient,  C  .   These 

'   m 

results  were  then  plotted  and  compared  with  the  experimental 
results  in  order  to  determine  which  value  of  C   gave  the 
best  fit.   These  curves  representing  the  "best  fit"  were 
then  plotted  on  Figure  (45)  and  were  compared  with  the  ex- 
perimental plots.   As  might  be  expected,  it  was  possible 
to  obtain  good  agreement  between  the  approximate  theory 
and  the  experimental  data  for  the  largest  depth  ratios.   In 
these  cases,  the  proximity  of  the  free  surface,  which  was 
disregarded  in  the  approximate  analysis,  was  probably  not 
particularly  important.   However,  for  the  shallowest  depth 
where  the  free  surface  was  more  influential,  it  appeared 
to  be  impossible  to  find  a  value  for  C  which  gave  reason- 
able agreement  with  the  plotted  experimental  results. 
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One  point  of  interest  in  Figure  (45)  was  to  compare 
the  two  curves  associated  with  a  maximum  depth  of  h/a=  6.0, 
the  one  case  corresponding  to  the  cylinder  located  on  the 
bottom,  d/a=  5.0,  and  the  other  with  the  cylinder  elevated 
one  radii  off  the  bottom,  d/a=  4.0.   As  can  be  seen  from 
the  figure,  the  latter  curve  is  considerably  below  the  for- 
mer and  corresponds  to  an  added  mass  coefficient  of  approxi 
mately  C  =  1.2.   The  data  corresponding  to  a  similar 
situation  with  the  cylinder  located  on  the  bottom  appears 
to  agree  relatively  well  with  the  curve  based  on  C   =  2.8. 
This  marked  increase  in  C   corresponding  to  the  case  where 
the  cylinder  was  located  on  the  bottom  apparently  results 
from  the  proximity  of  the  bottom. 

A  second  interesting  feature  in  Figure  (45)  was  the 
apparent  discontinuity  in  the  two  curves  associated  with 
the  two  shallowest  steps  near  2Tra/L^  0.3.   The  cause  of 
this  discontinuity  is  not  certain,  and  it  might  be  dis- 
missed as  scatter  in  the  experimental  data.   There  is, 
however,  a  second   explanation  which  is  equally  probable. 
It  is  possible  to  have  a  small  amount  of  reflection  at 
this  value  of  2Tra/L,  and  it  is  possible  that  this  reflec- 
ted waves  affected  the  wave  height  measurement.   Depending 
on  the  phase  of  the  reflected  wave,  it  may  either  add  or 
subtract  and  cause  an  apparent  increase  or  decrease  in  the 
value  of  the  wave  force  coefficient.   The  absence  of  this 
discontinuity  in  the  curves  corresponding  to  the  greater 
depths  is  in  agreement  with  this  hypothesis  since  reflec- 
tion decreases  with  relative  depth. 
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V.   CONCLUSIONS 

An  experimental  study  of  gravity  wave  interaction  with 
a  submerged  horizontal  cylinder  has  been  conducted.   From 
the  results  of  this  study,  the  following  conclusions  are 
drawn. 

1.  For  the  range  of  wave. heights  tested  in  the  experi- 
ment, the  horizontal  forces  resulting  from  gravity  wave 
interaction  with  the  cylinder  increased  linearly  with  wave 
height . 

2.  The  vertical  force  was  always  much  smaller  than 
the  horizontal  force  when  the  cylinder  was  located  on  the 
bottom.   The  downward  part  of  this  force  fluctuation  was 

crease  with  wave  amplitude.   At  very  small  amplitudes 
where  the  velocity  squared  contribution  to  the  pressure 
was  small,  a  sinusoidal  variation  of  the  vertical  force 
resulted. 

3.  For  the  case  corresponding  to  the  cylinder  sus- 
pended one  radii  off  the  bottom,  the  horizontal  force  was 
reduced  and  the  downward  force  was  greater  than  the  upward 
force . 

4.  For  cases  where  the  depth  of  submergence  was  large, 
the  approximate  analysis  appeared  to  give  valid  results 
provided  the  proper  value  of  the  added  mass  coefficient  was 
known . 
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5.   When  the  value  of  2-rra/L  becomes  small,  correspond- 
ing to  long  waves,  the  wave  force  trace  shows  features  dis- 
similar to  those  associated  with  the  more  sinusoidal  waves. 
In  particular,  a  null  period  in  the  horizontal  wave  trace 
occurs  during  the  passage  of  the  trough. 
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